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Abstract 

CO 

Using data collected with the LHCb detector in proton-proton collisions at a centre-of- 
mass energy of 7 TeV, the semileptonic decays B s — > DfXn~V and B s — > D°K + Xfi~V are 
detected. Two structures are observed in the D°K + mass spectrum at masses consistent 
with the known D sl (2536) + and -D* 2 (2573) + mesons. The measured branching fractions 
^ relative to the total B s semileptonic rate are B(B S — > D* S 2 X n~V) / 'B(B s — > Xji~V) = 

(3.3 ± 1.0 ± 0.4)%, and B(B° S -> D+X fx~V) / B(B° S ->■ X/i~]7) = (5.4 ± 1.2 ± 0.5)%, where 
the first uncertainty is statistical and the second is systematic. This is the first observation 
of the D* S 2 state in B s decays; we also measure its mass and width. 
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1 Introduction 



Much less is known experimentally about semileptonic B s decays, than for the lighter B 
mesons. In the case of the B s when the b — > c transition results in a single charm hadron 
this can be a D+, a D* + or another excited cs state. The relative proportion of these final 
states provides essential information on the structure of these semileptonic decays, and 
can be compared with QCD-based theoretical models. In this Letter we present a search 
for B s semileptonic decays, that might occur via an excited cs meson that disintegrates 
into final states containing D°K + . One such state is the Df ± , thought to be J p = 1 + , 
that decays into D*K, and another is the D*£, a possible 2 + state that has been observed 
to decay directly into DK [TJ. 

The LHCb detector [2] is a forward spectrometer constructed primarily to measure 
CP- violating and rare decays of hadrons containing b and c quarks. The detector elements 
are placed along the beam line of the LHC starting with the Vertex Locator (VELO), a 
silicon strip device that surrounds the proton-proton interaction region and is positioned 
8 mm from the beam during collisions. The VELO precisely determines the locations 
of primary pp interaction vertices, the locations of decays of long lived hadrons, and 
contributes to the measurement of track momenta. Other detectors used to measure 
track momenta comprise a large area silicon strip detector (TT) located before a 3.7 Tm 
dipole magnet, and a combination of silicon strip detectors (IT) and straw drift chambers 
(OT) placed afterward. Two Ring Imaging Cherenkov (RICH) detectors are used to 
identify charged hadrons. Further downstream an Electromagnetic Calorimeter (ECAL) 
is used for photon detection and electron identification, followed by a Hadron Calorimeter 
(HCAL), and a system consisting of alternating layers of iron and chambers (MWPC and 
triple-GEM) that distinguishes muons from hadrons (MUON). The ECAL, MUON, and 
HCAL provide the capability of first-level hardware triggering. 

In this analysis we use a data sample of approximately 20 pb" 1 collected from 7 TeV 
centre-of-mass energy pp collisions at the LHC during 2010. For the first 3 pb" 1 of these 
data a trigger was used that requires a single muon without any requirement that it 
misses the primary vertex, a trigger which was not available for the remainder of the data 
taking. This sample is well suited to determine the number of semileptonic B s decays, 
that we take as the sum of D+Xn~V, D°K + Xfi~V and D + K°Xfi~V decays, ignoring the 

small ~1% contribution from charmless B s decays. The entire 20 pb -1 sample, however, 
is useful for establishing signal significance, resonance parameter determination, and the 
ratio of numbers of events in the D°K + states. 

2 Selection criteria 

In both data samples backgrounds increase markedly with increasing track numbers. 
Thus, events are accepted only if the number of reconstructed tracks using the VELO 
and either the IT or OT is less than 100. Tracks were accepted based on similar criteria 
to those described in Ref. [2]. This results in only a 5.6% loss of signal in the 3 pb" 1 , and 
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a larger 9.4% loss over the entire 20 pb _1 sample. 

In this analysis we select a charm hadron that forms a vertex with an identified muon. 
We consider two cases: (i) Df — > K + K~ir + , that has a branching fraction of (5.50±0.27)% 

[1] - these are used to normalize the B s yield; (ii) D° — > K~ir + decays with a branching 
fraction of (3.89±0.05)% pp - these are combined with an additional K + that forms a 
vertex with the D° and the fi~ in order to search for B s semileptonic decays that might 
occur via an excited cs meson that decays into D°K + . In this Letter the mention of a 
specific final state will refer also to its charge-conjugate state. The selection techniques 
are similar to those used in a previous analysis [3J. Most charm hadrons are produced 
directly via pp — > ccX interactions at the LHC, where X indicates the sum over all other 
possible final state particles. We denote these particular charm reactions as "Prompt". 
Charm is also produced in pp — > bbX collisions where the fe-flavoured hadron decays into 
charm. These are called charm from b hadrons or "Dfb" for short. Muon candidates are 
selected using their penetration through the iron of the muon system. The candidates 
used in the analysis of the first 3 pb -1 sample must be those that triggered the event and 
have momentum transverse to the beam direction, p^, greater than 1200 MeV (we use 
units with c=l). 

The selection criteria for and D° mesons include identifying kaon and pion candi- 
dates using the RICH system. Cherenkov photon angles with respect to the track direction 
are examined and a likelihood formed for each particle hypothesis [2]. We also require 
that the pi of the kaons and pion be greater than 300 MeV, and that their scalar sum 
be greater than 2100 MeV (Df) or greater than 1400 MeV {D°). Since charm mesons 
travel before decaying, the kaon and pion tracks when followed backwards will most often 
not point to the primary vertex. The impact parameter (IP) is the minimum distance of 
approach of the track with respect to the primary vertex. We require that the x 2 formed 
by using the hypothesis that the IP is equal to zero, Xip, be > 9 for each track. The kaon 
and pion candidate tracks must also be consistent with coming from a common origin, the 
charm decay vertex, with vertex fit \ 2 P er number of degrees of freedom (ndof) < 6. This 
charm candidate's decay vertex must be detached from the closest primary interaction 
point. To implement this flight distance significance test we form a variable, Xfs' based 
on the hypothesis that the flight distance between the primary and charm vertices is zero, 
and require xls > Im- 
partial B s candidates formed from muon candidates must form a vertex with 
X 2 /ndof < 6, and point at the primary vertex: the cosine of the angle of the b pseudo- 
direction formed from the and muon vector momentum sum with respect to the 
line between the Dffi~ vertex and the primary vertex (cos 5) must be > 0.999. They 
must also have an invariant mass in the range 3.10 GeV< m(Dffi~) < 5.10 GeV. All of 
these requirements were decided upon by comparing the sidebands of the invariant mass 
distributions, representative of the background, with signal Monte Carlo simulation using 
PYTHIA 6.4 (5 event generation, and the GEANT4 [5] based LHCb detector simulation. 

The analysis for the D^X[i~V mode follows the same procedure as our previous 
D°Xfi~V study [3J, and uses the 3 pb _1 sample. The K + K~ir + mass spectra for both 
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the right-sign (RS K + K~tt + + and wrong-sign (WS K + K~n + + candidates, 
as well as the ln(IP/mm) distributions for events with mass combinations within ±20 
MeV of the D£ mass are shown in Fig. [1] for the pseudorapidity interval 2 < rj < 6. 
Here IP refers to the impact parameter of the Df candidate with respect to the primary 
vertex in units of mm. For both the RS and WS cases, we perform unbinned extended 
maximum likelihood fits to the two-dimensional distributions in K + K~n + invariant mass 
and ln(IP/mm), over a region extending from 80 MeV below the mass peak to 96 
MeV above. This fitting procedure allows us to determine directly the background shape 
from false combinations under the signal mass peak. The parameters of the Prompt 
IP distribution are found by examining directly produced charm [3J. The Monte Carlo 
simulated shape is used for the Dfb component. The fit separates contributions from Dfb, 
Prompt, and false combinations. The Prompt contribution is small. Background compo- 
nents for D* + — > 7T + D° — > ir + K + K~ and the reflection from A+ — > pK^ir + decay, where 
either a proton or a pion is wrongly identified as a kaon by the particle identification 
system, are also included. The shape of the D* + background is constrained to be equal 
to that of the — > K + K~n + signal peak and the yield is allowed to float, while the 
shape of the A+ reflection is determined from Monte Carlo and the yield is allowed to 
float within the uncertainty of our expectation. 

To evaluate more carefully the yield the fits are performed in t] bins and the 
detection efficiency in each bin is determined separately so as to remove uncertainty from 
differences in the t] dependent production observed in data compared to the Monte Carlo 
simulation. This procedure yields 2233±60 RS Dfb events in the D+X[i~V channel in 
the b pseudorapidity range 2 < 77 < 6, uncorrected for efficiency; the average detection 
efficiency is (1.07±0.03)%. This yield is then reduced by 5.1% for additional correlated b 
decay backgrounds as determined by simulation. 

3 Measurement of D°K + X fi~T7 

Semileptonic decays of B° s mesons usually result in a Df meson in the final state. It is 
possible, however, that the semileptonic decay goes to a cs excitation, which can decay 
into either DK or D*K resonances, or produces non-resonant DK. To search for these 
final states, we measure the D K + Xfi~u yield. To seek events with a D° candidate and 
an additional K + we require that the K + candidate has pt > 300 MeV, be identified as 
such in the RICH system, has X\p > 9? an d that the vector sum p^< of the D° and kaon 
be > 1500 MeV. The resulting partial B candidate must have an invariant mass in the 
range 3.09 GeV< m(D°K + < 5.09 GeV, form a vertex (x 2 /ndof < 3) and point at the 
primary vertex (cos<5 >0.999). In addition, we explicitly check that if the kaon candidate 
is assigned the pion mass and combined with the D°, it does not form a D* + candidate, 
by requiring the difference in masses m(A'~7r + 7r + ) — m(K~7r + ) — m(7r + ) > 20 MeV, in 
addition to the ±20 MeV requirement around the D° mass for m(K~ir + ). 

Figure Ufa) shows the D°K + invariant mass spectrum in the 3 pb _1 sample. D° 
candidates are chosen from K~7i + Xfi~V events with a K~ir + invariant mass within ±20 
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Figure 1: The invariant K + K~ir + mass spectra for events associated with a muon for 
the 3 pb _1 sample in the pseudorapidity interval 2 < t] < 6 for RS combinations (a) 
and WS combinations (c). Also shown is the natural logarithm of the IP distributions 
of the -D+ candidates for (b) RS and (d) WS muon candidate combinations. The 
labelling of the curves is the same on all four sub- figures. In descending order in (a): 
green-solid curve shows the total, the blue-dashed curve the Dfb signal, the black-dotted 
curve the sideband background, the purple-dot-dashed the misinterpreted — > pK~ir + 
contribution, the black dash-dash-dot curve the D* + — > tt + D° — > K + K~tt + contribution, 
and the barely visible red-solid curves the Prompt yield. The Dfb signal, the A+ reflection 
and D* + signal are too small to be seen in the WS distributions. The insert in (b) shows 
an expanded view of the region populated by Prompt charm production. 
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Figure 2: The mass difference m(K~n + K + ) — m(K~ir + ) added to the known D° mass 
for events with K~ir + invariant masses within ±20 MeV of the D° mass (black points) 
in semileptonic decays. The histogram shows wrong-sign events with an additional K~ 
instead of a K + . The curves are described in the text, (a) For the 3 pb _1 data sample 
and (b) for the 20 pb" 1 sample. 
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MeV of the D° mass. A clear narrow signal near threshold is seen corresponding to 
the _D s i(2536) + , but at a lower mass of 2392 MeV. An axial- vector state cannot decay 
into two pseudoscalar mesons but this resonance can decay into D*°K + . Since we do not 
reconstruct the 7 or 7r° from the D*°, the mass peak will be shifted down from its nominal 
value. However, because the resonance is so close to threshold, the mass resolution will 
still be very good resulting in a narrow peak. This final state was seen previously in B s 
semileptonic decays by the DO collaboration using — > D* + K§ decays [5]. There also 
appears to be a feature near the known mass of the _D* 2 (2573) + meson. The width of this 
state is not well measured; the PDG quotes 20±5 MeV pQ. Clearly there is a large excess 
over the wrong-sign background here evaluated using D°K~ mass combinations. 

In order to ascertain the size of the putative signals above background we perform an 
unbinned maximum likelihood fit. The data are fit with a threshold background function 
proportional to A4 p e~ aM , with A4 = m(D°K + ) — mo, where mo, the threshold point, is 
fixed at 2358.52 MeV. The fit determines p and a. We assume that the Tf > s — >• D° K + X fi~u 
signal above the background function is saturated by the and D* S 2 states. For the 
signal function we use a bifurcated Gaussian shape, whose relative widths above 
and below the peak are fixed from simulation. The mass and average width are fixed 
to the values 2391.6 MeV and 4.0 MeV, respectively, found using the higher statistics 
sample discussed below, while the simulation, including the effects of the missing D*° 
decay product, predicts a mass of 2392. 2±0. 3 MeV. The width is essentially due to the 
missing 7 or 7r° from the D*° to D° decay. There are 24.4±5.5 events. A relativistic 
Breit-Wigner signal shape convolved with the experimental resolution of 3.3 MeV (r.m.s.) 
is used in the region of the D*2 where both the mass and width are allowed to float in 
the fit. We find a mass value of 2559±9 MeV, a width of 24.1±9.2 MeV and 22.1±7.5 
events, where all of these uncertainties are statistical only. 

To confirm the _D*^~ signal we use the full data sample of 20 pb _1 , in which we accept all 

events that were triggered. While this sample is useful to increase statistics it suffers from 

a larger number of interactions per crossing, and multiple triggers, that makes it more 

— 

difficult to ascertain the total number of B s decays. The measurement of the relative 
yields of -D*^ to D^, however, will not be affected. Figure MJo) shows the resulting 
D°K + invariant mass spectrum. The difference between RS and WS events outside of 
the resonant peaks is consistent with background from other b decays as demonstrated 
by Monte Carlo simulation. We use the same fitting functions as above, but here we 
allow the mass and average width values of the bifurcated Gaussian to float while still 
fixing the ratio of widths above and below the peak from simulation. The fit to the 
yields 155±15 signal events, a D°K + mass of 2391.6±0.5 MeV, and 4.0±0.4 MeV for the 
width. For the _D*^~ we again allow the mass, the width and the number of events to float 
in the fit. We find a mass of 2569.4±1.6 MeV, a width of 12.1±4.5 MeV, and 82±17 
events. These errors are purely statistical. The previously measured mass and width 
values from the PDG are 2572.6±0.9 MeV and 20±5 MeV pQ. The probability of the 
background fluctuating to form the D*£ signal corresponds to eight standard deviations, 
as determined by the change in twice the natural logarithm of the likelihood of the fit 
without including this resonance and accounting for the change in the number of degrees 
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of freedom. 

The systematic uncertainty on the _D*^~ mass is determined from several calibration 
channels. For example, our measured D° mass differs from the known value by 0.2 MeV, 
though the known value has a 0.14 MeV error. We also see a variation on the order of 
0.3 MeV by varying the fit region and background shape, where we use a linear function 
instead of the threshold function. Thus we take ±0.5 MeV as the systematic uncertainty. 
We use the same method of changing the fits to find the systematic uncertainty on the 
width. The maximum observed change is 1.4 MeV. There is also a contribution from our 
uncertainty on the experimental resolution of ±0.5 MeV that contributes an additional 
0.7 MeV error on the width. Taking these two components in quadrature gives a width 
uncertainty of 1.6 MeV. 

The relative branching fractions are determined from the 20 pb _1 sample, assuming 
that the decays only into D*K final states, the _D*^~ decays only into DK final states, 
and isospin is conserved in their decays. Note that the only observed decays _D*^~ are 
to DK final states, while decays to D*K, although possible, have not yet been seen, 
including the study by the DO collaboration [5J. The D*£ / D^ event ratio is computed, 
correcting for the lower detection efficiency for Z)*^" of (0.516±0.017)%, compared with 
the D+ efficiency of (0.598±0.025)% as 

B(B° q -> D*tXu~V) 

- — - — - = 0.61 ±0.14 ±0.05. (1) 

B(B S D+Xn-V) 

The relative branching fraction of the D^ x with respect to the total B s semileptonic 
rate is measured using 24.4 ± 5.5 events in the 3 pb" 1 sample. The number of B s semilep- 
tonic decay events in this sample is evaluated from the efficiency corrected sum of the 
B s — > DfXfi~V events and twice the efficiency corrected B s — > D°XK + fi~17 yield. The 
efficiencies are 1.07% and 0.57%, respectively. The doubling of the D°K + X[i~V yield 
accounts for the missing D + K°X[i~V contribution, which is equal due to isospin symme- 
try. A small component of B — > D^KX[i~V is subtracted based on a branching fraction 
measurement from BaBar of (6.1 ±1.2) x 10~ 4 [7], reducing the D^XyTV yield by 3.2%. 
The overall uncertainty on the B s semileptonic yield is 6.6%. The main contributions 
to this error are the uncertainty on the absolute Df branching ratio of 4.9%, and the 
uncertainty on the amount of D°K + X fi~V events to add to the B s yield of 3.0%. The 
corresponding number for the D*£ branching fraction is computed also using this sample 
and the result from Eq. [TJ Correcting for the unreconstructed D + K° decays results in 
the doubling of the rates of the relative branching fractions, that we determine to be 



B(B° a -> DljX^V) 

B(lf 8 -> Xfi-V) 
B0i -> D+Xfx-V) 



-- (3.3 ±1.0 ±0.4)% 

= (5.4 ±1.2 ±0.5)%, (2) 
B(B" S -> X/i-v) 

where the systematic uncertainty for both includes a 5% error on the detection efficiency, 
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and the above mentioned 6.6% uncertainty on the number of B s semileptonic decays. In 
addition there is a systematic uncertainty of 8% on the D*£ yield estimated by varying 
the fit region, and background shape. Our branching fraction for the relative rate of 
decay is consistent with, but smaller than, the value of (9.8±3.0)% measured by DO [5J. 

4 Conclusions 

The first observation has been made of the rare semileptonic decay B° s — > 
D* 2 (2573) + X \x~V and its branching fraction relative to the total semileptonic B s de- 
cay rate has been measured as B(B S — > D*J^Xirv)j{B s — > X[i~V) = (3.3 ± 1.0 ± 0.4)%. 
For B s — > D s i(2536) + X [i~T> semileptonic decays the ratio is B(B S — >• Df 1 Xft~V)/B(B s — >• 
XyTV) = (5.4±1.2±0.4)%, where in both cases the first uncertainty is statistical and the 
second is systematic. We have assumed that the Df x decays only into D*K final states, 
the D*2 decays only into DK final states, and isospin is conserved in their decays. These 
values were predicted in the ISGW2 model as 3.2% and 5.7%, for D* S 2 and Df 1: respec- 
tively, in good agreement with our observations [8]. Another set of predictions based on 
the quark model are 1.8% and 2%, respectively [9]. The mass of the D*£ is measured to be 
2569.4±1.6±0.5 MeV, and the width as 12.1±4.5±1.6 MeV, in agreement with previous 
observations. 
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